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We present a development of a geometric approach to entanglement indicators. The method is 
applied to detect genuine multiqubit entanglement. The criteria are given in form of non-linear 
conditions imposed on correlation tensors. Thus they involve directly observable quantities, and 
in some cases require only few specific measurements to find multiqubit entanglement. The non- 
linearity of each of the criteria allows detection of entanglement in wide classes of states. In contrast 
to entanglement witnesses, which in the space of Hermitian operators define a hyperplane, the new 
conditions define a geometric figure encapsulating the non-fully entangled states within it. 
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I. INTRODUCTION 



Since quantum entanglement is both a basic resource 
in quantum communication and quantum information 
processing, and a fundamental phenomenon in consid- 
erations related to foundations of quantum physics, its 
qualitative and quantitative characterization becomes of 
a great importance for practical as well as purely theo- 
retical reasons 0, ■ 



II. REPRESENTATION OF STATES IN TERMS 
OF CORRELATIONS 



Any n-qubit state can be expressed as 
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where <7 Mfc G {1, a±, a%, 03} are the Pauli matrices of the 



nth observer. The coefficients T u 



are real numbers 



in [—1,1] given by correlation function values for mea- 
surements of products of Pauli operators 
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In contrast to the bipartite case, where the structure 
of entanglement is very simple (the state is either en- 
tangled or separable), for many subsystems the charac- 
terization of entanglement becomes more complex due 
to many possible ways of partitioning the whole system 
into subsystems. Several different approaches to detect 
genuine multipartite entanglement have been proposed: 
based on Bell inequalities [3-0|, using entanglement wit- 
nesses [1, HI , based on relations between elements of den- 
sity matrices [13, UU , utilizing Fisher information [12, HH 
and finally using correlation tensors [i~4| . 



In this contribution we further develop geometric ap- 
proach to entanglement detection proposed in Ref. [la ]. 
We show how this approach leads to necessary and suf- 
ficient conditions for various forms of multipartite en- 
tanglement and derive explicit criteria for different types 
of non-separability. The resulting conditions are of the 
form of non-linear combinations of correlation functions 
and sometimes the criteria are very simple and require 
only limited number of specific directly measurable data. 
The non-linearity of our criteria makes them often more 
versatile than entanglement witnesses. We provide exam- 
ples in which our criterion detects genuine multipartite 
entanglement of different families of quantum states. In 
our investigations, an inspiring role played the work of 
Yu et al. [16| , in which a general non-linear condition for 
two qubit entanglement was found. 



Tp u = (cr Ml ® ... <g) cr^Jp = Tr (p (T Ml <g) ... ® (T Mn ) . 

(2) 

A specific role is played by the components which involve 
only indices 1,2,3 (such indices will be denoted by Latin 
letters). The quantity: 
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where {e Im } 3 m=1 is a basis in K 3 , transforms like a ten- 
sor under local unitary transformations on the qubits. 
As a consequence, we refer to it as a correlation tensor 
of the state p. The whole object T = r Wr , i((> , where 
the indices take on values fXk = 0,1,2,3, will be called 
an extended correlation tensor. Its components with k 
zeros are (n — fc)-rank tensors. Extended correlation ten- 
sors belong to a real vector space with a natural scalar 
product 



(X,Y)=Y, X P Y P> 



(4) 



where (1 = (u 1; fj, n ) and uj, = 0, 1, 2, 3. We can gener- 
alize the notion of a scalar product between correlation 
tensors introducing a positive semidefinite metric G. A 
generalized scalar product has then the following form: 



(5) 
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This scalar product induces a G-norm: 

\\T\\ 2 G = (T,T) G . (6) 
Such a generalized scalar product and G-norm were used 




III. MULTIPARTITE ENTANGLEMENT 

Let us begin with a classification of entanglement of 
multipartite states. 

A pure n-partite state is called fc-product, if it can 
be represented as a tensor product of fc pure i m -partite 
states: 

\i>k- P ro&) = \i>ii)® ■■■®\il>i h ), (7) 

where of course, X)m=i * m = n - There can be different 
types of fc-product states corresponding to different ways 
of partitioning n into a sum of fc integers. We will refer 
to a definite type of fc-product state as {i\ + . . . + ik)- 
partition product state. Note that different fc-product 
states of the same type may involve in their partitions 
different physical subsystems. For example, (2 + 1)- 
partition product states of three particles A, B and G 
are: \ip A ) \iPbc), \^b) \*Pac) and \^ c ) \^ab)- We will re- 
fer to all this types of states as: (A + BC), (B + AC) 
and (AB + G) - product states respectively. 

The n-partite state p is called fc-separable, if it can 
be expressed as a probabilistic (convex) mixture of pure 
fc-product states: 

Pfc-scp = ^Pil^ft-prodKV'fc-prodl' ( 8 ) 

i 

In this terminology fully separable states are n-separable. 
If a state is not fc-separable, then it must involve entan- 
glement between at least n — k + 2 parties. Accordingly, 
a state is called genuinely multipartite entangled, if it is 
not biseparable, i.e. not 2-separable. 

Clearly, a set of fc-separable states is convex and it is 
a subset of (fc — Inseparable states as illustrated in Fig. 
[TJ They also touch each other in the following meaning: 
infinitesimally close to fully separable states there are 
states with entanglement between arbitrary number of 
subsystems. For an illustrative example see Section IVBl 



where T k ~ sep is a correlation tensor of some fc-separable 
state, and G is a metric in the sense of equation (JSJ) - It 
forms a sufficient condition for multipartite entanglement 
between at least n — k + 2 parties. 

The advantage of using scalar products to entangle- 
ment detection is that optimization over separable states 
in @ can be replaced with optimization over pure fc- 
product states only: 

max (T fc " sep , T) G = max ( V Pl Tf\ prod , T ) 

\ i / G 

< max (T k - prod ,T) G , (10) 

fk-prod 

This follows directly from linearity of a scalar product 
and convexity of fc-separable states. Since fc-separable 
states may involve different types of fc-product states, 
one should optimize over all possible partitions compat- 
ible with fc-separability. This reveals another feature of 
our condition of potential practical value: for different 
partitions n one can use different metrics to reveal that 
the state is not fc-product. All above implies the follow- 
ing modified condition: 

(V^c. max (T^,T) Gx < \\T\\ G ) 

\ rpk-prod J 

T is not fc-separable. (11) 

Finally, it is easy to show by adapting the reasoning of 
[l5j . that in case of any non-fc-separable multiqubit state, 
for any given partition, one can find such a metric G that 
LHS of the inequality in condition (jTTJ) holds, which leads 
to a necessary and sufficient condition for genuine multi- 
partite entanglement. Indeed, the condition of rejecting 
full separability, originally formulated in (l5| in terms of 
density matrices, can be put as follows: 

max Tr (pGp pro d) < Tr(pGp) =>■ p is not product state, 

Pprod 

(12) 

where now G is a positive semidefinite superoperator. In 
15] it is shown, that, conversly, if a state p is not sepa- 
rable, there exists a positive semidefinite superoperator 
G, such that the inequality (TT2l holds. To state this fact 



IV. NECESSARY AND SUFFICIENT 
CONDITIONS FOR MULTIPARTITE 
ENTANGLEMENT 

To indicate a case of a non-fc-separability we shall use 
the following simple geometrical observation, which is a 
corollary of the results of [l5j |. namely that: 



3 G max(T fc -^,T) G < 

rpk-sep 



\T\ 



G 



T is not fc-sep. 



(9) 




FIG. 1. Sets of fc-separable states: all sets are convex, the 
darker is the color of a set the more separable are its states, 
i.e. a set of fc-separable states contains as its subsets all more- 
than-fc-separable states. 
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clearly let us introduce new notation. Let {/ M }^ =0 de- 
note standard basis in the space of 2 x 2 matrices over 
complex numbers M2(C), that is the space of operators 
acting on a Hilbert space C 2 of a single qubit. Any den- 
sity matrix p can then be decomposed as p = ^2pppfp, 
where p = {pi, . . . , p„} and fp, = / M1 ® . . . ® / Mn is a ba- 
sis in a tensor product space (M^)®" = M2". According 
to [HI, if a state p is not separable, then there exists a 
superoperator G, such that: 



maxTr(pGp prod ) < Tr(pGp). 



(13) 



In the case of rejecting fc-separability instead of full sep- 
arability we proceed in full analogy, except for the fact 
that for optimilizing the LHS of (IT21 for each partition, 
we treat all elements of this partition as a single system 
(of course it can be of higher dimension) . Having k sin- 
gle systems we can find a metric G for each partition 
separately, since correctness of condition (fTB")) does not 
depend on dimensions of single systems under consider- 
ation. By choosing a specific basis (putting = er^), 
all above considerations can be translated into correla- 
tion tensor representation of a state ((T|). Indeed, having 
a density matrix of any state in correlation tensor form 
, we can find its components in a standard basis in the 
following manner: 



PP 



E 



(14) 



The scalar product between density matrices p and p' in 
metric G can now be reformulated as a scalar product 
between extended correlation tensors: 



Tr(pGp') = (p,p')c 



7 t 7 j T^Up^GppTgUpg 

7<5 



(15) 



Above transformation implies, that the metric operator 
G which defines a scalar product of correlation tensors 
(cf. ([5])) corresponding to a scalar product of density 
matrices with superoperator G has the following form: 



(16) 



Altogether we have shown the following proposition 
fully characterizing ^-separability: 

Proposition 1 An n-particle state endowed with ex- 
tended correlation tensor T is not k-separable if and only 
if for every partition n into k subsystems there exists a 
metric G^, such that the following inequality holds: 



max (T^^T)^ < (T,T) G „ 



(17) 



V. EXAMPLES 

We present here a few examples of applications of 
Proposition [1] It leads to simple sufficient conditions for 
multipartite entanglement that are sometimes also neces- 
sary and detects entanglement of various classes of states 
(which is impossible using entanglement witnesses). 



A. Three-qubit entanglement 

To reveal a genuine 3-partite entanglement in a three 
qubit state we have to exclude the case of biseparability. 
We derive several sufficient criteria using different metric 
tensors. 



1. Standard metric 

Here we shall give a condition which is unbiased in 
its formulation with respect to any family of entangled 
states. It uses the diagonal metric, that with Gpp in the 
form of a Kronecker delta Sp?. 

Let T be a correlation tensor of a 3-qubit state, and 
let T 2+1 be a correlation tensor of a (2 + l)-partition 
product 3-qubit state. Assuming standard (Euclidean) 
scalar product in space K 9 : 



3 

E ' x ' 



A?: 7 h ij' 



ijk , 



(18) 



we show, that max„ maxj,2+i (T 2+1 ,T), where a denotes 
permutation of biproduct states (here three possible 
splittings), is upper bounded by: 



max max . 



E (I^K 11 *) _ T cr(22i) \ + \T a (33i)\) j 
i=l 

(19) 

where o~{0 ® O', 1) means, that the second maximiza- 
tion is done over local orthogonal transformations ap- 
plied to subsystems over which, for a given a, summa- 
tion is not performed, e.g. if a is a trivial permutation, 
CT(0(g)6',i) = 6® 6' ®L 

Clearly f 2+1 = f 2 ®^. Any pure two qubit state can 
be expressed in a Schmidt basis in the form: cos (9 1 00) + 
sin#|ll), which has the following nonvanishing terms of 
correlation tensor T 2 : 



Tn = sin 20 
T 22 = -sin26> 
T 33 = l 

while T 1 is a Bloch vector: 



fx 



[mi,m 2 ,m 3 j, 



with 



(20) 



(21) 



4 



Hence the scalar product (T 2+1 ,T), where a denotes 
some permutation of indices defining to which subsys- 
tems tensors T 2 and T 1 correspond, can be expressed 
as: 

3 

(f*+\f) = J2 {(Ta(iu) - T a[22i) ) sin 20 + T ff(334) ) to,. 

i=l 

(22) 

Since to, is arbitrary unit vector, the maximization over 
m gives: 

max(f (T 2 + 1 ,f) = 



(22i)) Sin 26* + ^(334)) 

\ i=l 



(23) 



Finally we have to maximize over 8, over permutations a 
of subsystems and over all possible local rotations O®0' 
applied to subsystems, over which we do not sum in Eq. 
(f2"3")l . At this stage we do not need to maximize over 
rotations applied to the subsystem over which we sum in 
the above equation, since maximization over a unit to is 
equivalent to maximization over all possible rotations of 
this vector, and the maximization over rotations of the 
scalar product (T 2+1 ,T) can be performed over any of 
the tensors in this product. 

Adopting the notation introduced in Eq. (|19p and us- 
ing the fact that for reals T\ and r<i we have (r\ sin 28 + 



r 2 Y< n 



1 7*2 1) j we reach the following condition: 



Proposition 2 If the following inequality holds: 



<r,a(O0O' 



^ (\ T cr(lli) - T a (22i) | + 1^(334)1) < 11^1 



(24) 

then the state described by correlation tensor T is gen- 
uinely 3-partite entangled. 

However, as the left-hand side of this proposition may be 
strictly bigger than Eq. (|2"3")l , in some cases it is more 
effective to directly use condition 



maxmax(T 2+1 ,T) < ||T| 

a ^2+1 



(25) 



For example, for a generalized 3-partite GHZ state: 

\GHZ a ) = cosa|000) + sina|lll), (26) 
mixed with white noise: 



1 



v\GHZ a )(GHZ a \ +(l-v)-l 



(27) 



Eq. (|23p is maximized for vanishing 8, and reads 

3 sin 2a. The same value is obtained after ap- 
plying local rotations as we verified numerically. Since 



the squared length of the correlation tensor equals i> 2 (T 
3 sin 2 2a), we find that if v exceeds the critical value 



VI + 3 sin 2 2a ' 



(28) 



the state is genuinely 3-partite entangled for any a. Lo- 
cal rotations applied to the left-hand side of (|24p make its 
value higher than the one which follows from Eq. (|23[) . 
and accordingly Proposition [2] does not detect as many 
states as direct application of Eq. (|2"3"|) . Finally, note 
that the critical visibility (|28l) holds for all possible lo- 
cally unitarily equivalent 3-partite GHZ states. For a 
symmetric GHZ state (that is for a = tt/4), v cr u = \. 
Note that even for an arbitrarily small but finite a, the 
state (|2"6"|) is 3-partite entangled, while for a — it is 
fully separable. 

In case of a 3-partite W state \W 3 ) = -^=(|100) + |010) + 

1 001)) mixed with white noise: 

^(p) = »|W3><W 3 | + (l-«)^l (29) 

the critical visibility for detection of genuine 3-qubit en- 
tanglement with the condition (|2~4"|) has been calculated 
numerically, and is equal to v cr i t ~ 0.636. 

We stress that this approach is quite versatile, it al- 
lows to detect three-particle entanglement despite the 
fact that GHZ and W states are of a different nature. 



2. Modified metric and generalized Schmidt decomposition 
of the correlation tensor 

The criterion in Proposition [2] can be made more effi- 
cient and simplified by changing a metric and applying a 
generalized Schmidt decomposition [l7| to the correlation 
tensor. 

In order to arrive at a simple criterion, we modify the 
metric (hence also scalar product in the sense of Eq. ([SJ)) 
to get rid of terms of the T^^j) type. Thus, we put 



17-112 _ T 2 -V^T 

I I \ mod / j ijk / j c 

i.j.k 1 I 



2 

<r(33Z)' 



(30) 



Using this metric, the inequality (j2~4"]l can be rewritten 
as: 



max max 



i) \ i= i 



EC- 



(iii) 



er(22i) 



| 2 < ||f,,2 



mod ' 



(31) 

Since the condition Q is valid in any metric G, we have 
the following modified criterion: 

Proposition 3 If the following inequality holds: 



max max 



T 



)) <\\T 



|2 

I modi 



(32) 
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where \\f\\ 2 mod = J2lj,k=i T ijk ~ E?=i T a(33i)' then the 
state described by correlation tensor T is genuinely 3- 
partite entangled. 

We can further simplify this condition applying handy 
features of a generalized Schmidt decomposition to the 
correlation tensor [l7j. According to Theorem 1 in [l7j . 
for any tensor 



T = T, 



e" , with jfe = 1, ...,d, (33) 



which detects a smaller class of entangled states, but is 
experimentally very handy. One can measure compo- 
nents of T which enter ||T|| mod , and once the sum ((30|) 
is above 2, a genuine three particle entanglement is con- 
firmed. 

All above analysis can be performed in full analogy 
in case of 4-partite states. Due to complexity of formu- 
las, conditions for 4-qubit entanglement are presented in 
Appendix A. 



where {e^}f m=1 is a basis in some d-dimensional vector 
space, there exists a basis ® ... <g> sf n ,ik = 1, ...,d, 
which we shall call a generalized Schmidt basis, in which 
the components i of tensor T: 



have the following properties: 

•T^. w) =0forl<i<i<d, 

• ... i„i s non- negative if 
at most one of the if. differs from d, 



(34) 



(35) 



(36) 



Ti.J > K,_ in \ if J < i r for all r = 1 ".(37) 



B. GHZ metric 

In this section we shall study an approach that favours 
a certain family of states, namely the GHZ ones. 

Consider the so-called Greenberger-Horne-Zeilinger 
state of n qubits: 

\GHZ) = ^=(\0...0) + \l...l)). (39) 

It has the following non- vanishing elements of extended 
correlation tensor [18[ : 

Ty...y x .... x = ( — 1) , k = 0, 1, |_' 2i 2^J 



Assume that the correlation tensor T in (f3"2"| is expressed 
in a generalized Schmidt basis. Then the property ([33)) 
implies, that from among the following two groups of 
terms of the correlation tensor: 

{Tiiii T\\2, T113} 



{T221, T222, T223} 

only one term in each group is nonzero. Let us assume 
without loosing generality that Tin is the maximal gen- 
eralized Schmidt coefficient and only T221 is non-zero. 
This implies that for all a: 



(Hi) — T a (22i) 



< ITiuI + ITmiI <2|T m | 



This property gives the following criterion: 

Proposition 4 // the following inequality holds: 

I ITI I 2 ~v 9T 

IK II mod ^J-maxj 



(38) 



where \ \f\\ 2 mod = Tf jk - J2Li T *(33i)> and T " 



^2 

is the maximal possible value a correlation tensor element 
for the given state, then the state described by correlation 
tensor T is genuinely 3-partite entangled. 

The above condition can be simplified further to a weaker 
one: 



I ^11 mo< j > 2, 



r GHZ 

-Z...Z0...0 



(40) 



and with indices permuted. In turns out that Proposition 
[T]with a diagonal metric G, 



(X,Y) G =J2XpGpY p , 



(41) 



for the specific case of Gp = IT™ 2 ! and additionally 
putting Go,...,o = 0, leads to optimal detection of a gen- 
uine multipartite entanglement of noisy GHZ states 



P = v\GHZ)(GHZ\ + (l-v) — t. 



(42) 



Using the metric, the right-hand side of Proposition Q] 
for state p is given by (2" — l)v 2 . In order to find the 
maximum of the left-hand side over bi-product states, we 
write 



_ ^rp rpbi— prod* 



a = v 



GHZ Nrrnbi— prod 



peguz 



where we introduce a convenient notation for summing 
over non-zero elements of a "GHZ" metric, and note that 
after transforming correlation tensors to density opera- 
tors this reads 

L/v = 2"- 1 (p 1 , 1 + p 1<2 n + p 2 „ ! + p 2 »,2«) - 1 (44) 
= 2"Tr(p GHZ p bi - prod ) - 1. (45) 

In the last step we use Tr(p GHZ p bi " prod ) < |, which holds 
for all bi-separable states [l9j ] . Therefore, L < (2™ -1 — l)v 
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and the state p is shown to be genuinely multipartite 
entangled if 



v > 



1 



2™ - 1 



(46) 



which is known to be optimal [111, ]20| . 

The same metric can reveal genuine multipartite en- 
tanglement of many other states. Consider generalized 
GHZ states 



\GHZ a ) =cosa|0...0) +sina |1 . . . 1) , 



(47) 



with a € [0, The non- vanishing components of its 
correlation tensor are given by permutations of indices of 
the following ones 



r, 



y--v- 



T, 



Z...Z0...0 



(-l) fc sin2a, fe = 0, 1, 



1 for k even, 

cos 2a for k odd. 



(48) 



Taking again the "GHZ" metric, one can repeat the proof 
which led to (|46|) with the only difference that now 



Tr(p GH5 V i - prod ) <cos 2 a. 

In this way we obtain that generalized GHZ state mixed 
with white noise is genuinely multipartite entangled for 



v > 



2 n cos 2 a - 1 



1 



(49) 



Finally, note that this state is fully separable only for 
a = 0. Already for infinitcsimally small a it can involve 
entanglement between all n parties. Clearly, a similar 
statement would hold for a generalized GHZ state be- 
tween n — 1 parties and the nth party having an uncor- 
related state. Therefore, in infinitesimal neighborhood 
of the state |0 . . . 0) there are states with entanglement 
between an arbitrary number of subsystems. 



VI. CONCLUSIONS 

We have presented several sufficient criteria for a multi- 
partite entanglement in the form of nonlinear conditions 
imposed on correlations of the tested state. The condi- 
tions are given in a convenient and simple form, and can 
be directly applied to given families of entangled states. 

An important advantage of our criteria is that in many 
cases only few definite measurements suffice to detect 
multiqubit entanglement. 

Presented criteria are more general than entanglement 
witnesses due to their nonlincarity. A single new cri- 
terion detects a genuine entanglement of many different 
families of states, whereas one definite witness can detect 
entanglement of one family of states only. 

Here we have given only several examples, however one 
can construct infinitely many other ones. Note, that only 



our conditions with the GHZ metric involved correlations 
of all qubits as well as only some of them. Note that this 
is the case for the universal two qubit entanglement con- 
dition given in (l6| . thus this seems to be a promising 
direction of a further research. Different series of condi- 
tions of such a kind, with surprising properties, will be 
presented elsewhere [211 ] . 
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Appendix A 

1. Genuine 4- partite entanglement in four-qubit 

states 

a. Exclusion of biseparability 

To exclude a biseparability of a 4-qubit state, we have 
to verify the condition ^ for the case of maximizing over 
(3 + 1)- and (2 + 2)-partition product states: 



max(T 3+1 ,T) G < 



ITI 



and ( max(T 2+2 ,T) G 

V T 2 + 2 
rp rpbisep 



< \\T\ 



(Al) 



where G denotes a metric operator in a vector space in 
which 4-qubit correlation tensors are embedded. Using 
a version of a " GHZ" metric for a four-qubit system, in 
which only T .( 112 2)-type terms occur: 



I 1 \\GHZ = J llll t ^1122 t J 1221 

I rri'2 _t rj^2 I /T~i2 

t J 1212 ~r ^2121 + J 2112 



rpl 

' ± T1\\ 

rpl 

J 2222? 



(A2) 



we calculate the first term in (IA1|) . Since T 3+1 is a pure 
state, we have f 3+1 = f 3 ® f 1 . Taking f 1 = m = 



mi,m 2 



m 3 ], with \J m\ + i 



1, and assuming 
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T 3 ® T 1 is (ABC + D)-type product, one obtains: 

(f 3 <g> f\ f) = TmnnTun + Tii 2 m 2 r 1122 
+ Ii 22 miTi 22 i + r 22 imiT 22 n 
+ ri 2 im 2 Ti 2 i 2 + T 2 i 2 miT 2 i 2 i 
+ 7 2 iim 2 T 2 n 2 + T 222 m 2 T 2222 . (A3) 

Due to properties (|33|) . (l3f)|) and (|37p of a generalized 
Schmidt decomposition [1 7j j applied now to quantum 
states, any 3-qubit pure state can be expressed as: 

\ip) = cos(wi)|000) +cos(w 2 )sin(a;i)|001} 
+ e i0 cos(w 3 )sin(a;i)sin(w 2 )|OlO) 
+ cos(w4) sin(wi) sin(a; 2 ) sin(w3)|100) 
+ sin(wi) sin(w 2 ) sin(a;3) sin(w4)|Hl)- (A4) 



In this parametrization, terms of T 3 occuring in (|A3|) 
have the following form: 

T\\2 = T\2\ = Tan = r 222 = 
T\\\ = — Ti 22 = — T 22 i = — T 2 i 2 

= sin(2wi) sin(w 2 ) sin(cj 3 ) sin(w 4 ). (A5) 

Hence (|A3|) simplihcs to: 

(T 3 <g f l ,f) — (T lm - T 1221 - T 22U - T 2121 ) 

xrni sin(2wi) sin(w 2 ) sin(w3) sin(w4). 

(A6) 



The maximization of (|A6|) over mi, Wi, cj 2 , W3, W4 is triv- 
ial. Let us denote local orthogonal transformations 
0\ (g O2 O3 (g O4 as Otot . We finally obtain: 

max(T 3+1 ,T) = max|T m i - Xi 22 i - T 22 n - T 2i2 i| . 

^ 3+1 6 tot 

(A7) 

In complete analogy we can find inequalities for other 
types of (3 + l)-partition product states (ABD+C, 
ACD+B, A+BCD), which leads to the following inequal- 
ities: 

max iTim — Ti 22 i — T 22 n — T2121 1 < H^Hghz 

o tot 

max |Tim — T 22 n — Ti 2 i 2 — T 2 n 2 | < HrHcHz 

Otot 

max |T lm — T 1122 — T 2121 — T 2112 | < II^IIghz 

Otot 

max |Tim — Tn 22 — Ti 22 i — Ti 2 i 2 | < ||T||q HZ . 

Otot 

(A8) 

In case of pure states these inequalities allow us to check 
if given state is (3+l)-partition product or not: 



Proposition 5 // all the inequalities \A8\) hold, then the 
pure state described by correlation tensor T is not (3 + 1)- 
partition product. 



Now we have to calculate the second element of the 
conjunction in (| Al|) involving maximization over (2+2)- 
partition product states. Since 



max(f 2+2 , f) = max (T 2 <g f ' 2 , f ), (A9) 

T 2 + 2 t 2 T' 2 



we need to explicitely express T 2 (g T 2 . This is very 
simple due to equations (|20|) : 



(f 2 : 


in 


= -sin26> 


(f 2 : 


122 


= sin 29 


(t 2 : 


)33 


= 1 


(T 2 : 


111 


= - sin(26»' 


(f 2 : 


)22 


= sin(20') 


if' 2 ; 


)33 


= 1. 



(A10) 



The only nonvanishing terms of the tensor T 2 (g T 2 are: 



{f 2 <i 


5T 2 )mi 


= sin26>sin26>' 


(f 2 a 


5 T 2 )n 22 


= - sin 29 sin 29' 


(f 2 a 


5 T 2 )n33 


= -sin26> 


(f 2 a 


5 T 2 ) 22 n 


= -sin26»sin26»' 


(f 2 a 


5 T 2 ) 2222 


= sin26»sin26» / 


(f 2 a 


3 T 2 ) 22 33 


= sin 29 


(f 2 a 


5 T 2 )33H 


= - sin 29' 


(f 2 a 


5 T 2 )33 22 


= sin 29' 


(f 2 a 


5 T 2 )3333 


= 1. 



(All) 



Substituting these terms one obtains (in case of maxi- 
mizing over (AB+CD)-product states): 

max (f 2 (gf' 2 ,f) = 

max (sin 29 sin 29' (Tim - Tn 22 - T 22 n + T 2222 ) 

6tot,B,9' 



sm29(T 2233 -T 1133 ) 

sin 29' (T33 22 — T3311) + 13333)- 



(A12) 



Since only 4 terms of (|A1 1|) occur in GHZ metric, the 
expression (|A12[) has, for maximizing over (AB+CD)- 
type product states, the following simplified form: 

max (f 2 (gf' 2 ,f) = 

rp2 rp'2 

= max (sin20sin(20')(J 1 iiii-Tii 22 -T 22 ii+T 2222 )) 

Otot-.O.O' 

= max |Tim - Tn 22 - T 2211 + T 2222 |. (A13) 

Otot 

Taking into account other types of (2+2)-partition prod- 
uct states (that is of type (AC+BD) and (AD+BC)) we 
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obtain the following set of inequalities: 



max 

Otot 


Tim 


— -M.122 


— 12211 


+ T2222 


< 


max 

Otot 


Tim 


— 7l2i2 


— ?2121 


+ ?2222 1 


< 


max 

6 tot 


Tim 


— T.221 


~ T2112 


+ T2222 1 


< 



(A14) 



Finally, we obtain the Proposition, which is a direct con- 
sequence of condition (|Al[) : 



Proposition 6 If all the inequalities HA8\) and \A 14\ ) 

hold, then the state described by correlation tensor T is 
genuinely 4-partite entangled. 



b. Exclusion of 3 -separability 

Since there is only one type of 4-partite 3-product 
state, that is (2 + 1 + l)-partition product, the condi- 
tion ([9|) has the following form: 



max(T 2+1+1 ,T) G < \\f\\ 2 G 



T2+1 + 1 



f^f 3sep . (A15) 



be expressed as: 

(T 2+1+1 ,T) = 
3 

= E ((^(ihj) _T '<T(22ij))sin26' + T (T (3 3ii ))m i rt l . 

(A17) 

Now we use Cauchy-Schwartz inequality: 

3 



E {(Ta(iiij) - T CT (2 2 i,))sin20 + T a ( 33ij) ) m,n, 



< 



tr(33ij)j 



'■,i=i 



E^ 2 ^ 

3 

E {(T a (iuj) -7V(22ij))sin20-|-T CT (33i J ))' 



2/3 



< 



^=1 



E (p* 

lly) - T a{22ij)) sin 2 # + ^(33^)) | ( Al 8) 

The last inequality follows from Cauchy-Schwartz in- 
equality and the fact that: 



We proceed analogously to the case of excluding bisep- 
arability of 3-partite state: T 2+1+1 = f 2 <g> f 1 ® T" 1 , 
and we choose Schmidt basis for T 2 , in which the only 
nonvanishing terms are: 



Tn = sin26» 
T 22 = - sin 29 
T33 = 1, 



(A16) 



while T 1 = m = [mi, 7712,7713], with ^/ m 2 + m 2 , + m| = 
1, and T' 1 = n = [711,712,713], with v' rl i + n 2 + n i = 1- 
Hence the scalar product (T 2+1+1 ,T), where cr denotes 
proper permutation of indices refering to subsystems, can 



Em 2 <l^E™^ L 



From now on we can proceed directly as in the case of 
maximizing over (2 + l)-partition product states obtain- 
ing: 

Proposition 7 If the following inequality holds: 



maxmax ~ T a(22i 3 ) \ + \ T a(33ij) I) 2 < \\T \ 

° tot \ij=i 

(A19) 

then the state described by correlation tensor T is bisep- 
arable or genuinely multiqubit entangled. 
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